Measurement of the Negative Muon Anomalous Magnetic Moment to 0.7 ppm 
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The anomalous magnetic moment of the negative muon has been measured to a precision of 
0.7 parts per million (ppm) at the Brookhaven Alternating Gradient Synchrotron. This result is 
based on data collected in 2001, and is over an order of magnitude more precise than the previous 
measurement of the negative muon. The result a^^- = 11 659 214(8)(3) x 10^^" (0.7ppm), where the 
first uncertainty is statistical and the second is systematic, is consistent with previous measurements 
of the anomaly for the positive and negative muon. The average for the muon anomaly is a;i(exp) = 
11659 208(6) X lO"^'^ (0.5 ppm). 

PACS numbers: 13.40.Em, 12.15.Lk, 14.60.Ef 



The anomalous magnetic moments of the muon and 
the electron have played an important role in the develop- 
ment of the standard model. Compared to the electron, 
the muon anomaly has a relative sensitivity to heavier 
mass scales which typically is proportional to (m^/mg)^. 
At the present level of accuracy the muon anomaly gives 
an experimental sensitivity to virtual W and Z gauge 
bosons as well as a potential sensitivity to other, as yet 
unobserved, particles in the few hundred GeV/c^ mass 
range l]. 

We report our result for the negative muon anomalous 
magnetic moment a^- = {g — 2) /2 from data collected in 
early 2001. The measurement is based on muon spin pre- 
cession in a magnetic storage ring with electrostatic fo- 
cusing. The same experimental technique was used as in 
our most recent measurements of and a similar 

precision of 0.7 ppm was achieved. Detailed descriptions 
of the apparatus may be found elsewhere 0, 0, IE 0. • 

For polarized muons moving in a uniform magnetic 
field B perpendicular to the muon spin and to the plane 
of the orbit and in an electric quadrupole field E, which 
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is used for vertical focusing the angular frequency 
difference, oja between the spin precession frequency and 
the cyclotron frequency, is given by 
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The dependence of uja on the electric field is eliminated by 
storing muons with the "magic" 7 = 29.3 9], which cor- 
responds to a muon momentum p = 3.09 GeV/c. Hence 
measurement of LOa and of B, in terms of the free proton 
NMR frequency ujp and the ratio of muon to proton mag- 
netic moments A, determines a^. At the magic 7, the 
muon lifetime is approximately 64.4 fis and the (g — 2) 
precession period is 4.37 /iS. With a field of 1.45 T in our 
storage ring Q, the central orbit radius is 7.11 m. 

The difference frequency uJa was determined by count- 
ing the number N{t) of decay electrons above an energy 
threshold. The time spectrum of decay electrons is then 
given by 



N{t) = iVoe"*/(^") [1 + Asm{iUat + 0a)] 



(2) 



The normalization A'o, asymmetry A, and phase (j)a vary 
with the chosen energy threshold. 

The measurement of the magnetic field frequency LUp is 
based on proton NMR in water. A trolley with 17 NMR 
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Normal Skew 

Quad -0.28 0.11 
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FIG. 1: A two-dimensional multipole expansion of the 2001 
field averaged over azimuth from one out of 20 trolley mea- 
surements. Half ppm contours with respect to a central az- 
imuthal average field Bo = 1.451 269 T are shown. The multi- 
pole amplitudes relative to Bo are given at the beam aperture, 
which had a radius of 4.5 cm and is indicated by the circle. 



TABLE L Systematic uncertainties for the Up analysis. 



Source of uncertainty 


Size [ppm] 


Absolute calibration of standard probe 


0.05 


Calibration of trolley probe 


0.09 


Trolley measurements of Bq 


0.05 


Interpolation with fixed probes 


0.07 


Uncertainty from muon distribution 


0.03 


Others'^ 


0.10 


Total systematic error on Up 


0.17 



^ higher multipoles, trolley temperature and voltage re- 
sponse, eddy currents from the kickers, and time-varying 
stray fields. 



probes was moved typically every three days throughout 
the entire muon storage region. About 150 fixed NMR 
probes distributed around the ring in the top and bottom 
walls of the vacuum chamber were used to interpolate the 
field between trolley measurements. The system was cali- 
brated with respect to a standard probe with a spherical 
H2O sample The homogeneity of the field in 2001 
(Figure ^1 was similar to that achieved for the opposite 
polarity field in 2000 0. 

The field (B) weighted with the analyzed event sam- 
ple was obtained from two largely independent anal- 
yses, whose results were found to agree to within 
0.05 ppm. Its final value is expressed in terms of 
the free proton resonance frequency and is given by 
u)p/{2t:) = 61 791 400(11) Hz (0.2 ppm). Table U lists the 
uncertainties. The improved 2001 uncertainties resulted 
from refinements in the calibration measurements, and 
from an upgraded system to determine the azimuthal 
trolley position in the storage ring. 

The 2001 LUa data taking was similar to that in 2000. 
However, the hardware energy threshold of the detectors 
was kept lower and equal for all counters at 0.9 GeV com- 
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FIG. 2: The Fourier spectrum of the residuals of a fit to the 
five free parameters given in Eq.|21for the high (top) and low 
(bottom) n-value data. The corresponding CBO frequencies, 
located at 491 kHz (top), and 419 kHz (bottom) as well as 
their {g — 2) sidebands are clearly visible. Dashed lines indi- 
cate the (g — 2) frequency. 



pared to 1.0-1.4 GeV in 2000. This was made possible 
by reducing the intensity of the injected beamjwhich in 
turn reduced the light flash in the detectors 0, |^ . These 
factors allowed all the detectors to be turned on and be 
stable by 32 /is after beam injection, as opposed to 50 fj,s 
in 2000. As a result of the reduced rates, the fraction 
of overlapping signals (pilcup) after 32 /is in 2001 was 
comparable to the pileup fraction after 50 /is in 2000. In 
2000 the field focusing index n, which is proportional to 
the electric field gradient, was n = 0.137, corresponding 
to a horizontal coherent betatron oscillation frequency 
(CBO) of 466 kHz Q. This frequency was close to twice 
the {g — 2) frequency of 229 kHz, which resulted in a 
sizable uncertainty in the fitted uJa value Q. In 2001 
we used two different ri-values, n = 0.122 and n — 0.142, 
which resulted in CBO frequencies, 419 kHz and 491 kHz 
that are further from twice the ((7 — 2) frequency (see Fig- 
ure [5J. Consequently, the uncertainty caused by CBO is 
smaller. Furthermore it also reduced the correlation be- 
tween the CBO and detector gain effects in the fits to the 
time spectrum. 

Two independent implementations of the algorithm 
to reconstruct the electron times and energies from the 
calorimeter signals were used. The frequency iOa was de- 
termined by fitting the time distribution of decay elec- 
trons. Five independent analyses were performed in or- 
der to probe the systematic uncertainties and, of course. 
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to protect against mistakes. All five results agreed within 
the expected statistical deviations due to different data 
selection and weightings. These analyses are described 
below. 

Two of the analyses used slightly different parametriza- 
tions that included CBO modulations and fitted the 
combined electron spectrum in the energy range 1.8-3.4 
GeV. In the third analysis, the counts were weighted with 
the experimentally determined energy-dependent modu- 
lation asymmetry, which optimized the statistical power 
of the data. This method permitted the analyzed en- 
ergy range to be extended. We used an energy range 
of 1.5 to 3.4 GeV, which together with the asymmetry 
weighting resulted in a 10% improvement of the statisti- 
cal uncertainty. As in the first two analyses, the resulting 
spectrum of weighted counts was fitted to a function that 
parametrized all known and statistically significant per- 
turbations. 

The remaining analyses fit the ratio 0, Q formed by 
randomly assigning the data to four statistically inde- 
pendent subsets ni to n^. The subsets were rejoined 
in u{t) = ni{t) + n2{t) and v{t) = n^{t — Ta/2) + 
ni{t + Ta/2), where Ta is an estimate of the {g — 2) 
period, and then combined to form the time spectrum 
r{t) = [u{t) - v{t)\l\u{t) + v(t)\. The {g - 2) rate modu- 
lation of V is 180° degrees out of phase compared to that 
of u, and to sufficient precision can be described by 
A'smiiOot + (\)a)- The ratio r(t) is largely insensitive to 
changes of observed counts on time scales larger than 

Ta = 2TT/uJa ~ 4/XS. 

In one of the ratio analyses, the sensitivity to CBO was 
reduced by combining the data from both n-values and 
all detectors prior to fitting. The data were fitted from 
32 yus after injection when all detectors were on. In the 
second ratio analysis the data were fitted separately for 
each calorimeter and n— value. The fits began between 
24 jjLS and 32 /is, and required the parametrization of the 
CBO effects in the fit function. 

Changes in the radial and vertical muon distributions 
with time were quantified, and were found to have neg- 
ligible effect on LDa- A small reduction in the pulsed 
electrostatic quadrupole voltages ^] during the measure- 
ment period could change the vertical muon distribution. 
Analysis of the data from scintillator counter hodoscopes 
placed in front of the calorimeters combined with a beam 
tracking calculation and a GEANT based simulation set 
a systematic error limit of 0.03 ppm. The muon radial 
distribution is determined by the magnetic field and the 
momentum distribution [2, Ui| . The magnetic field does 
not change with time after injection, except due to the 
field from eddy currents induced by the fast kicker 0- 
This was measured, and found to have a negligible effect 
on the muon radial distribution. Muons of lower mo- 
menta decay earlier in the laboratory frame than muons 
of higher momenta. The momentum distribution of the 
stored beam thus changes during the 600 /iS measure- 
ment period. The effect on uia due to this change was 
studied in simulation, and was found to be 0.03 ppm. 




FIG. 3: Comparison of the uia values from the five analyses 
for the low-n (filled) and high-n (open) data sets. Analysis 4 
used only the combined low and high-n data (square). The 
divisions on the vertical axis are separated by 1 ppm, and 
the indicated uncertainties are statistical. The systematic 
uncertainties are considerably smaller. 



TABLE II: Systematic uncertainties for the combined Ua anal- 
ysis. 



Source of errors 


Size [ppm] 


Coherent betatron oscillations 


0.07 


Pileup 


0.08 


Gain changes 


0.12 


Lost muons 


0.09 


Others'!" 


0.11 


Total systematic error on Ua 


0.21 



^ AGS background, timing shifts, E field and vertical os- 
cillations, beam debunching/randomization, binning and 
fitting procedure. 



The results for uja for the two n-values are consistent, 
see Figure 131 and were combined for each of the analyses. 
The values for uja from the five analyses are in agreement 
to within variations expected from the differences in the 
analyzed event samples and the treatment of the data. 
The analysis techniques are expected to have somewhat 
different sensitivities to different systematic effects. De- 
tailed comparisons of the results, using all analyzed data 
as well as only the data in overlap, showed no evidence 
for unaccounted systematic differences. The five result- 
ing values for oja were combined in a simple arithmetic 
mean to obtain a single value for uja- 

The resulting frequency value is uja/{2T:) = 
229 073.59(15) (5) Hz (0.7ppm), which includes a correc- 
tion of -f 0.77(6) ppm for contributions to Eq.^ caused 
by vertical oscillations (0.30 ppm) and for the effect of 
the horizontal electric fields on muons with 7 7^ 29.3 
(0.47 ppm). The stated uncertainties account for strong 
correlations among the individual results, both statistical 
and systematic. Table HTl lists the systematic uncertain- 
ties in the combined result with these correlations taken 
into account. 

After the LUp and cua analyses were finalized separately 
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and independently, was evaluated. 



The result is 



R 



X-R 



11659 214(8)(3) x 10"^" (0.7ppm) 



(3) 

where i?^- = uia/^p = 0.003 707 208 3(2 6) and A = 
= 3.183 345 39(10) (Hi]. This new result 
is in good agreement with the average of = 
0.003 707 204 8(2 5) as predicted by the CPT theorem. 



The difference AR = R„ 



R,, 



(3.5 ± 3.4) X 10- 



The new average is i?^ = 0.003 707 206 3(2 0) and 

a^(exp) = 11659 208(6) X 10^1° (0.5ppm), (4) 

in which the total uncertainty consists of 5 x 10^^° 
(0.4 ppm) statistical uncertainty and 4 x 10~^° (0.3 ppm) 
systematic uncertainty. The correlation of systematic un- 
certainties between the data sets has been taken into ac- 
count. The combined result for the positive muon 
a^+(exp) = 11659 203(8) x 10"!° (0.7 ppm) has a statis- 
tical uncertainty of 6 x 10"^" (0.6 ppm) and a systematic 
uncertainty of 5 x 10^^" (0.4 ppm). It is shown in Fig- 
ure ^together with the new result for the negative muon 
and their average. 
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FIG. 4: Measurements of by E821 with the SM predic- 
tions (see text for discussion) . Uncertainties indicated on the 
measurements are total uncertainties. 

The standard model prediction for consists of QED, 
hadronic and weak contributions. The uncertainty on the 
standard model value is dominated by the uncertainty 
on the lowest-order hadronic vacuum polarization. This 
contribution can be determined directly from the anni- 
hilation of e'^e~ to hadrons through a dispersion inte- 
gral [l^ . The indirect determination using data from 
hadronic t decays, the conserved vector current hypoth- 
esis, plus the appropriate isospin corrections, could in 
principle improve the precision of ap(had). However, 
discrepancies between the r and the e'^e~ results ex- 
ist [TsL H^ . The two data sets do not give consis- 
tent results for the pion form factor. Using e~^e~ an- 
nihilation data the corresponding theoretical value is 
a^(SM) = 11659181(8) x lO-^" (0.7ppm). The value 
deduced from r decay is larger by 15 x 10^^" and has a 
stated uncertainty of 7 x 10""'^'' (0.6 ppm). The differ- 
ence between the experimental determination of and 
the standard model theory using the e+e^ or r data for 
the calculation of the hadronic vacuum polarization is 
2.7 a and 1.4 a, respectively. 

This is the final analysis of the anomalous magnetic 
moment from experiment E821 at the Brookhaven Alter- 
nating Gradient Synchrotron. We aim to substantially 
improve our result in a new measurement and look for- 
ward to continued efforts to improve the theoretical eval- 
uation. 
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